Retirado do artigo Miller et al. (2019). Distance sampling in R. Journal of Statistical Sofware 89(1)
# instalar pacotes necessários
#install.packages("Distance")
# instalar pacotes adicionais
#install.packages("mrds")
#install.packages("dsm")
#install.packages("mads")
#install.packages("dsims")
# carregar pacotes
library(Distance)
library(dplyr)
library(DT)
library(flextable)
library(ggplot2)
library(lubridate)
library(plotly)
library(readr)
library(readxl)
library(stringr)
library(tibble)
library(tidyr)
# carregar as funções da pasta R
# carregar função script_carregar_funções_pasta_r.R
source(
paste0(
here::here(),
"/R/minhas_funcoes.R"
)
)
# carregar dados
cutia_tap_arap <- transformar_para_distanceR_covariaveis() |>
filter(
Region.Label == "Resex Tapajos-Arapiuns",
sp_name == "Dasyprocta croconota"
) |>
drop_na(distance)
# readr::write_excel_csv(
# cutia_tap_arap,
# paste0(
# here::here(),
# "/data/cutia_tap_arap.csv"
# ),
# )
cutia_tap_arap <- cutia_tap_arap |>
select(
Region.Label,
Area,
Sample.Label,
Effort,
distance
)
cutia_tap_arap |>
DT::datatable(filter = "top")
dados_completos |>
group_by(
uc_name,
ea_name
) |>
filter(day_effort == max(day_effort)) |>
drop_na(distance) |>
ungroup()
dados_completos |>
filter(
uc_name == "Resex Tapajos-Arapiuns",
sp_name == "Dasyprocta croconota"
) |>
datatable(filter = list(position = "top"))
Variáveis necessárias para o data.frame:
Region.Label: vetor fator com o estrato contendo o
transecto (pode ser uma estratificação pré-amostragem - UCs - ou
pós-amostragem - ex. região, estado, bioma)
Area: vetor numérico contendo a área do
estrato;
Sample.Label: vetor númerico contendo a identidade
(ID) do transecto
object: nome adicional, ver seção 6;
detected: nome adicional, ver seção 6;
Effort: vetor númerico contendo o esforço do
transecto (para linhas seu comprimento, para pontos o número de vezes
que o ponto foi visitado)
size: vetor numérico copntendo o tamanho do grupo
observado;
distance: vetor numérico de distâncias
observadas;
Month:
OBs:
Sp:
mas:
HAS:
Study.Area:
Transectos que foram amostrados, mas que não tiveram observações (n =
0) devem ser incluídos no conjunto de dados com NA nas
observações de distância e qualquer outra covariael para a qual não se
tenha observação.
# cutia_tap_arap |>
# complete(Region.Label, Sample.Label, sp_name) |>
# datatable(filter = list(position = "top"))
Jogar a imputacao de NAs pra dentro da funcao carregar
dados completos.
# desenha o grafico com a distribuicao de distancias perpendiculares
cutia_tap_arap |>
filter(distance >= 1,
distance <= 14) |>
plotar_distribuicao_distancia_interativo()
summary(cutia_tap_arap$distance)
Ajustando um modelo ao dados das cutias Dasyprocta
croconota, configurando uma distância limite de 20m e usando
Half-normal como key function usando o argumento
key, sem termo de ajuste.
cutia_tap_arap_filtrado <- cutia_tap_arap |>
filter(distance >= 1,
distance <= 14)
# ajustando a função de detecção para uma distancia de truncamento de 20, 15, 10 e 5 metros
# dsitancias de truncamento
dist_truncamento <- list(
`14 metros` = 14,
`12 metros` = 12,
`10 metros` = 10,
`7 metros` = 7
)
# Key function - Half-normal
cutia_tap_arap_hn <- purrr::map(
dist_truncamento,
\(.x) ajuste_modelos_distance_hn(
cutia_tap_arap_filtrado,
truncamento = .x
)
)
cutia_tap_arap_hn
Ajustando um modelo ao dados da cutia Dasyprocta croconota,
configurando uma distância limite de 20m e usando Hazard rate
como key function usando o argumento key.
# ajustando a função de detecção para uma distancia de truncamento de 20, 15, 10 e 5 metros
# Key function - Hazard-rate
cutia_tap_arap_hr <- purrr::map(
dist_truncamento,
\(.x) ajuste_modelos_distance_hr(
cutia_tap_arap_filtrado,
truncamento = .x
)
)
Ajustando um modelo ao dados das cutias Dasyprocta
croconota, configurando uma distância limite de 20m e usando
Uniform como key function usando o argumento
key, sem termo de ajuste.
# ajustando a função de detecção para uma distancia de truncamento de 20, 15, 10 e 5 metros
# Key function - Uniform
cutia_tap_arap_unif <- purrr::map(
dist_truncamento,
\(.x) ajuste_modelos_distance_unif(
cutia_tap_arap,
truncamento = .x
)
)
summarize_ds_models(
cutia_tap_arap_hn$`14 metros`$`Sem termo`,
cutia_tap_arap_hn$`14 metros`$Cosseno,
cutia_tap_arap_hn$`14 metros`$`Hermite polinomial`,
cutia_tap_arap_hr$`14 metros`$`Sem termo`,
cutia_tap_arap_hr$`14 metros`$Cosseno,
cutia_tap_arap_hr$`14 metros`$`Polinomial simples`,
cutia_tap_arap_unif$`14 metros`$Cosseno,
cutia_tap_arap_unif$`14 metros`$`Polinomial simples`
)
summarize_ds_models(
cutia_tap_arap_hn$`12 metros`$`Sem termo`,
cutia_tap_arap_hn$`12 metros`$Cosseno,
cutia_tap_arap_hn$`12 metros`$`Hermite polinomial`,
cutia_tap_arap_hr$`12 metros`$`Sem termo`,
cutia_tap_arap_hr$`12 metros`$Cosseno,
cutia_tap_arap_hr$`12 metros`$`Polinomial simples`,
cutia_tap_arap_unif$`12 metros`$Cosseno,
cutia_tap_arap_unif$`12 metros`$`Polinomial simples`
)
summarize_ds_models(
cutia_tap_arap_hn$`10 metros`$`Sem termo`,
cutia_tap_arap_hn$`10 metros`$Cosseno,
cutia_tap_arap_hn$`10 metros`$`Hermite polinomial`,
cutia_tap_arap_hr$`10 metros`$`Sem termo`,
cutia_tap_arap_hr$`10 metros`$Cosseno,
cutia_tap_arap_hr$`10 metros`$`Polinomial simples`,
cutia_tap_arap_unif$`10 metros`$Cosseno,
cutia_tap_arap_unif$`10 metros`$`Polinomial simples`
)
summarize_ds_models(
cutia_tap_arap_hn$`7 metros`$`Sem termo`,
cutia_tap_arap_hn$`7 metros`$Cosseno,
cutia_tap_arap_hn$`7 metros`$`Hermite polinomial`,
cutia_tap_arap_hr$`7 metros`$`Sem termo`,
cutia_tap_arap_hr$`7 metros`$Cosseno,
cutia_tap_arap_hr$`7 metros`$`Polinomial simples`,
cutia_tap_arap_unif$`7 metros`$Cosseno,
cutia_tap_arap_unif$`7 metros`$`Polinomial simples`
)
O que tem que ter?
Os gráficos (probabilidade de detecção pela distância, com a curva ajustada, exemplo abaixo, fazer no ggplot), resultado do goodness of fit (gof_ds()), cada modelo vai ter que ter um nome diferente numa tabela(?)
plot(cutia_tap_arap_hn, breaks = seq(0, 20, 2.5))
plot(cutia_tap_arap_hn_herm, breaks = seq(0, 20, 2.5))
plot(cutia_tap_arap_hn_cos, breaks = seq(0, 20, 2.5))
plot(cutia_tap_arap_hr, breaks = seq(0, 20, 2.5))
plot(cutia_tap_arap_hr_poly, breaks = seq(0, 20, 2.5))
plot(cutia_tap_arap_hr_cos, breaks = seq(0, 20, 2.5))
Podemos usar a função summary para obter informações
importantes sobre o modelo.
lista_modelos <- list(
cutia_tap_arap_hn,
cutia_tap_arap_hn_herm,
cutia_tap_arap_hn_cos,
cutia_tap_arap_hr,
cutia_tap_arap_hr_poly,
cutia_tap_arap_hr_cos
)
purrr::map(lista_modelos, \(x) summary(x))
summarize_ds_models(
cutia_tap_arap_hn,
cutia_tap_arap_hn_herm,
cutia_tap_arap_hn_cos,
cutia_tap_arap_hr,
cutia_tap_arap_hr_poly,
cutia_tap_arap_hr_cos
)
O resultado inclui detalhes sobre o dado e a especificação do modelo, assim como dos coeficientes (\(\beta_{j}\)) e sua inceteza, a média do valor de detectabilidade e sua incerteza e uma estimativa da abundância na área coberta pela amostragem (sem levar em consideração o tamanho dos agrupamentos, ou bandos).
Para visualizar quão bem a função de detecção se ajusta aos dados quanto temos as distâncias exatas podemos usar um plot de quantis empíricos x teóricos (Q-Q plot). Ele compara a função de distribuição cumulativa (CDF) dos valores ajustados da função detecção a distribuição empírica dos dados (EDF).
Também podemos usar o teste de Cramér-von Mises para testar se os pontos da EDF e da CDF tem origem na mesma distribuição. O teste usa a soma de todas as distâncias entre um ponto e a linha y = x para formar a estatística a ser testada. Um resultado significativo fornece evidência contra a hiipótese nula, sugerindo que o modelo não se ajusta bem aos dados.
# ajustando um modelo Half-normal
cutia_hn <- ds(data = cutia_tap_arap_15,
truncation = 20,
transect = "line",
key = "hn",
adjustment = NULL)
# conduzindo o teste dfe bondadede ajuste de Cramer-von Mises
gof_ds(cutia_hn)
gof_ds(cutia_hr_time)
O resutlado do teste aponta que o modelo Half-normal deve ser descartado.
Testes de bondade de ajuste de chi-quadrado são gerados usando a
função gof_ds quando as distâncias forneceidas estão
categorizadas.
Uma vez que temos um conjunto de modelos plausíveis, podemos utilizar
o cirtériode informaçãode Akaike (AIC) para selecionar entre os modelos
o que melhor se ajusta aos dados utilizando a função
summarize_ds_models.
# gerando uma tabela de seleção de modelos usando AIC
summarize_ds_models(cutia_hn, cutia_hr_time, cutia_hr_time_size)
O melhor modelo é o Hazard-rate com tempo de senso e tamanho do grupo como covariáveis.
Para obter a abundância na região de estudo, primeiro calculamos a abundância na área amostrada para obter \(N_c\) e em seguida escalonamos esse valor para toda a área de estudo multiplicando \(N_c\) pela razão entre a área amostrada e a área da região. Para estimar a abundância na área amostrada, utilizamos as estimativas de probabilidade de detecção no estimador de Horvitz-Thompson.
Quando fornecemos os dados no formato correto (“flatfile”)
ds irá automaticamente calcular as estimativas de
abundância baseado nas informações de amostragem presenta nos dados.
summary(cutia_hn)
Summary statistics: fornece as áreas, aŕea de amostragem, esforço, número de observações, número de transectos, taxa de encontro, seus erros padrões e coeficientes de variação para cada estrato;
Abundance: fornece estimativas, erros padrões, coeficientesde variação, intervalos de confiança inferior e superior, graus de liberdade para a estimativa de abundância de cada estrato;
Densidade: lista as mesmas estatísticas de Abundance, só que para densidade.
contar_n_repeticoes_trilha() - conta o número de vezes
que cada trilha foi visitada
Ajuste Hermite pollynomial usa od código "herm"
e polinomial simples "poly".
Podemos incluir covariáveis utilizando o argumento
formula = ~ .... Abaixo, está especificado um modelo
“Hazard-rate” para os dados de cutia q ue inclui o tempo de senso como
covariável e uma distância limite de 20m.
cutia_hr_time <- cutia_tap_arap_15 |>
ds(truncation = 20,
key = "hr",
formula = ~ cense_time)
Adicionando uma segunda covariável: tamanho do grupo.
cutia_hr_time_size <- ds(data = cutia_tap_arap_15,
truncation = 20,
transect = "line",
key = "hr",
formula = ~ cense_time + size)
plot(cutia_hr_time)
plot(cutia_hr_time_size)
cutia_esec_terra_meio <- transformar_para_distanceR_covariaveis() |>
filter(
Region.Label == "Esec da Terra do Meio",
sp_name == "Dasyprocta croconota"
) |>
drop_na(distance)
# desenha o grafico com a distribuicao de distancias perpendiculares
cutia_esec_terra_meio |>
filter(distance >= 1,
distance < 15) |>
plotar_distribuicao_distancia_interativo(largura_caixa = 1)
Ajustando um modelo ao dados das cutias Dasyprocta
croconota, configurando uma distância limite de 20m e usando
Half-normal como key function usando o argumento
key, sem termo de ajuste.
cutia_esec_terra_meio_filtrado
cutia_esec_terra_meio_filtrado <- cutia_esec_terra_meio |>
filter(distance >= 1,
distance < 15)
# ajustando a função de detecção para uma distancia de truncamento de 20, 15, 10 e 5 metros
# dsitancias de truncamento
dist_truncamento <- list(
#`20 metros` = 20,
`15 metros` = 15,
`12 metros` = 12,
`10 metros` = 10
)
# Key function - Half-normal
cutia_esec_terra_meio_hn <- purrr::map(
dist_truncamento,
\(.x) ajuste_modelos_distance_hn(
cutia_esec_terra_meio_filtrado,
truncamento = .x
)
)
cutia_esec_terra_meio_hn
Ajustando um modelo ao dados da cutia Dasyprocta croconota,
configurando uma distância limite de 20m e usando Hazard rate
como key function usando o argumento key.
# ajustando a função de detecção para uma distancia de truncamento de 20, 15, 10 e 5 metros
# Key function - Hazard-rate
cutia_esec_terra_meio_hr <- purrr::map(
dist_truncamento,
\(.x) ajuste_modelos_distance_hr(
cutia_esec_terra_meio_filtrado,
truncamento = .x
)
)
Ajustando um modelo ao dados das cutias Dasyprocta
croconota, configurando uma distância limite de 20m e usando
Uniform como key function usando o argumento
key, sem termo de ajuste.
# ajustando a função de detecção para uma distancia de truncamento de 20, 15, 10 e 5 metros
# Key function - Uniform
cutia_esec_terra_meio_unif <- purrr::map(
dist_truncamento,
\(.x) ajuste_modelos_distance_unif(
cutia_tap_arap,
truncamento = .x
)
)
summarize_ds_models(
cutia_esec_terra_meio_hn$`20 metros`$`Sem termo`,
cutia_esec_terra_meio_hn$`20 metros`$Cosseno,
cutia_esec_terra_meio_hn$`20 metros`$`Hermite polinomial`,
cutia_esec_terra_meio_hr$`20 metros`$`Sem termo`,
cutia_esec_terra_meio_hr$`20 metros`$Cosseno,
cutia_esec_terra_meio_hr$`20 metros`$`Polinomial simples`,
cutia_esec_terra_meio_unif$`20 metros`$Cosseno,
cutia_esec_terra_meio_unif$`20 metros`$`Polinomial simples`
)
summarize_ds_models(
cutia_esec_terra_meio_hn$`15 metros`$`Sem termo`,
cutia_esec_terra_meio_hn$`15 metros`$Cosseno,
cutia_esec_terra_meio_hn$`15 metros`$`Hermite polinomial`,
cutia_esec_terra_meio_hr$`15 metros`$`Sem termo`,
cutia_esec_terra_meio_hr$`15 metros`$Cosseno,
cutia_esec_terra_meio_hr$`15 metros`$`Polinomial simples`,
cutia_esec_terra_meio_unif$`15 metros`$Cosseno,
cutia_esec_terra_meio_unif$`15 metros`$`Polinomial simples`
)
summarize_ds_models(
cutia_esec_terra_meio_hn$`10 metros`$`Sem termo`,
cutia_esec_terra_meio_hn$`10 metros`$Cosseno,
cutia_esec_terra_meio_hn$`10 metros`$`Hermite polinomial`,
cutia_esec_terra_meio_hr$`10 metros`$`Sem termo`,
cutia_esec_terra_meio_hr$`10 metros`$Cosseno,
cutia_esec_terra_meio_hr$`10 metros`$`Polinomial simples`,
cutia_esec_terra_meio_unif$`10 metros`$Cosseno,
cutia_esec_terra_meio_unif$`10 metros`$`Polinomial simples`
)
summarize_ds_models(
cutia_esec_terra_meio_hn$`12 metros`$`Sem termo`,
cutia_esec_terra_meio_hn$`12 metros`$Cosseno,
cutia_esec_terra_meio_hn$`12 metros`$`Hermite polinomial`,
cutia_esec_terra_meio_hr$`12 metros`$`Sem termo`,
cutia_esec_terra_meio_hr$`12 metros`$Cosseno,
cutia_esec_terra_meio_hr$`12 metros`$`Polinomial simples`,
cutia_esec_terra_meio_unif$`12 metros`$Cosseno,
cutia_esec_terra_meio_unif$`12 metros`$`Polinomial simples`
)
cutia_parna_serra_pardo <- transformar_para_distanceR_covariaveis() |>
filter(
Region.Label == "Parna da Serra do Pardo",
sp_name == "Dasyprocta croconota"
) |>
drop_na(distance)
# desenha o grafico com a distribuicao de distancias perpendiculares
cutia_parna_serra_pardo |>
filter(distance < 15,
distance > 0) |>
plotar_distribuicao_distancia_interativo()
Ajustando um modelo ao dados das cutias Dasyprocta
croconota, configurando uma distância limite de 20m e usando
Half-normal como key function usando o argumento
key, sem termo de ajuste.
# ajustando a função de detecção para uma distancia de truncamento de 20, 15, 10 e 5 metros
# dsitancias de truncamento
dist_truncamento <- list(
`20 metros` = 20,
`15 metros` = 15,
`10 metros` = 10,
`5 metros` = 5
)
# Key function - Half-normal
cutia_parna_serra_pardo_hn <- purrr::map(
dist_truncamento,
\(.x) ajuste_modelos_distance_hn(
cutia_parna_serra_pardo,
truncamento = .x
)
)
cutia_parna_serra_pardo_hn
Ajustando um modelo ao dados da cutia Dasyprocta croconota,
configurando uma distância limite de 20m e usando Hazard rate
como key function usando o argumento key.
# ajustando a função de detecção para uma distancia de truncamento de 20, 15, 10 e 5 metros
# Key function - Hazard-rate
cutia_parna_serra_pardo_hr <- purrr::map(
dist_truncamento,
\(.x) ajuste_modelos_distance_hr(
cutia_parna_serra_pardo,
truncamento = .x
)
)
Ajustando um modelo ao dados das cutias Dasyprocta
croconota, configurando uma distância limite de 20m e usando
Uniform como key function usando o argumento
key, sem termo de ajuste.
# ajustando a função de detecção para uma distancia de truncamento de 20, 15, 10 e 5 metros
# Key function - Uniform
cutia_parna_serra_pardo_unif <- purrr::map(
dist_truncamento,
\(.x) ajuste_modelos_distance_unif(
cutia_parna_serra_pardo,
truncamento = .x
)
)
summarize_ds_models(
cutia_parna_serra_pardo_hn$`20 metros`$`Sem termo`,
cutia_parna_serra_pardo_hn$`20 metros`$Cosseno,
cutia_parna_serra_pardo_hn$`20 metros`$`Hermite polinomial`,
cutia_parna_serra_pardo_hr$`20 metros`$`Sem termo`,
cutia_parna_serra_pardo_hr$`20 metros`$Cosseno,
cutia_parna_serra_pardo_hr$`20 metros`$`Polinomial simples`,
cutia_parna_serra_pardo_unif$`20 metros`$Cosseno,
cutia_parna_serra_pardo_unif$`20 metros`$`Polinomial simples`
)
summarize_ds_models(
cutia_parna_serra_pardo_hn$`15 metros`$`Sem termo`,
cutia_parna_serra_pardo_hn$`15 metros`$Cosseno,
cutia_parna_serra_pardo_hn$`15 metros`$`Hermite polinomial`,
cutia_parna_serra_pardo_hr$`15 metros`$`Sem termo`,
cutia_parna_serra_pardo_hr$`15 metros`$Cosseno,
cutia_parna_serra_pardo_hr$`15 metros`$`Polinomial simples`,
cutia_parna_serra_pardo_unif$`15 metros`$Cosseno,
cutia_parna_serra_pardo_unif$`15 metros`$`Polinomial simples`
)
summarize_ds_models(
cutia_parna_serra_pardo_hn$`10 metros`$`Sem termo`,
cutia_parna_serra_pardo_hn$`10 metros`$Cosseno,
cutia_parna_serra_pardo_hn$`10 metros`$`Hermite polinomial`,
cutia_parna_serra_pardo_hr$`10 metros`$`Sem termo`,
cutia_parna_serra_pardo_hr$`10 metros`$Cosseno,
cutia_parna_serra_pardo_hr$`10 metros`$`Polinomial simples`,
cutia_parna_serra_pardo_unif$`10 metros`$`Sem termo`,
cutia_parna_serra_pardo_unif$`10 metros`$Cosseno,
cutia_parna_serra_pardo_unif$`10 metros`$`Polinomial simples`
)
summarize_ds_models(
cutia_parna_serra_pardo_hn$`5 metros`$`Sem termo`,
cutia_parna_serra_pardo_hn$`5 metros`$Cosseno,
cutia_parna_serra_pardo_hn$`5 metros`$`Hermite polinomial`,
cutia_parna_serra_pardo_hr$`5 metros`$`Sem termo`,
cutia_parna_serra_pardo_hr$`5 metros`$Cosseno,
cutia_parna_serra_pardo_hr$`5 metros`$`Polinomial simples`,
cutia_parna_serra_pardo_unif$`5 metros`$`Sem termo`,
cutia_parna_serra_pardo_unif$`5 metros`$Cosseno,
cutia_parna_serra_pardo_unif$`5 metros`$`Polinomial simples`
)
purrr::map_df(
list(
cutia_esec_terra_meio_hn$`20 metros`,
cutia_esec_terra_meio_hr$`20 metros`
),
\(.x) purrr::map_df(.x, \(.y) summarize_ds_models(.y))
)
purrr::map_df(
cutia_esec_terra_meio_hn$`15 metros`,
\(.x) summarize_ds_models(.x)
)
purrr::map_df(
cutia_esec_terra_meio_hn$`10 metros`,
\(.x) summarize_ds_models(.x)
)
purrr::map_df(
cutia_esec_terra_meio_hn$`5 metros`,
\(.x) summarize_ds_models(.x)
)
Fonte: datuopinion.com
sagui_mont_tumuc <- transformar_para_distanceR_covariaveis() |>
filter(
Region.Label == "Parna Montanhas do Tumucumaque",
sp_name == "Saguinus midas"
) |>
drop_na(distance)
sagui_mont_tumuc |>
plotar_distribuicao_distancia_interativo()
Warning: Continuous y aesthetic
ℹ did you forget `aes(group = ...)`?
sagui_mont_tumuc_hn <- sagui_mont_tumuc |>
ajuste_modelos_distance_hn(truncamento = 10)
Fitting half-normal key function
AIC= 353.849
Starting AIC adjustment term selection.
Fitting half-normal key function
AIC= 353.849
Fitting half-normal key function with cosine(2) adjustments
AIC= 354.52
Half-normal key function selected.
Starting AIC adjustment term selection.
Fitting half-normal key function
AIC= 353.849
Fitting half-normal key function with Hermite(4) adjustments
AIC= 355.818
Half-normal key function selected.
sagui_mont_tumuc_hr <- sagui_mont_tumuc |>
ajuste_modelos_distance_hr(truncamento = 10)
Fitting hazard-rate key function
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).AIC= 355.417
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).Starting AIC adjustment term selection.
Fitting hazard-rate key function
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).AIC= 355.417
Fitting hazard-rate key function with cosine(2) adjustments
Warning: Detection function is not strictly monotonic!AIC= 356.358
Hazard-rate key function selected.
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).Starting AIC adjustment term selection.
Fitting hazard-rate key function
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).AIC= 355.417
Fitting hazard-rate key function with simple polynomial(4) adjustments
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).AIC= 360.38
Hazard-rate key function selected.
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).
sagui_mont_tumuc_hn |>
purrr::map(\(.x) plot(.x))
$`Sem termo`
NULL
$Cosseno
NULL
$`Hermite polinomial`
NULL
sagui_mont_tumuc_hr |>
purrr::map(\(.x) plot(.x))
$`Sem termo`
NULL
$Cosseno
NULL
$`Polinomial simples`
NULL
summarize_ds_models(
sagui_mont_tumuc_hn$`Sem termo`,
sagui_mont_tumuc_hn$Cosseno,
sagui_mont_tumuc_hn$`Hermite polinomial`,
sagui_mont_tumuc_hr$`Sem termo`,
sagui_mont_tumuc_hr$Cosseno,
sagui_mont_tumuc_hr$`Polinomial simples`
)
sagui_mont_tumuc_hn |>
purrr::map(\(.x) gof_ds(model = .x))
$`Sem termo`
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 0.301163 p-value = 0.134156
$Cosseno
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 0.301163 p-value = 0.134156
$`Hermite polinomial`
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 0.301163 p-value = 0.134156
sagui_mont_tumuc_hr |>
purrr::map(\(.x) gof_ds(model = .x))
$`Sem termo`
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 0.288212 p-value = 0.145969
$Cosseno
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 0.288212 p-value = 0.145969
$`Polinomial simples`
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 0.288212 p-value = 0.145969
# área de estudo, tamanho da área de estudo, area coberta pelo esforço amostral, esforço amostral em metros, número de detecções, número de transectos (ea), taxa de encontro, coeficiente de variação da taxa de encontro
sagui_mont_tumuc_hn$`Sem termo`$dht$individuals$summary[1:9]
# área de estudo, tamanho da área de estudo, trilhas ou estações amostrais, esforço total em cada trilha, abundância estimada em cada estação amostral, número de detecções em cada estação amostral, área total amostrada
sagui_mont_tumuc_hn$`Sem termo`$dht$individuals$Nhat.by.sample[1:8]
# total, densidade estimada, erro padrão da densidade destimada, coeficiente de variação da densidade destimada, intervalo de confiança inferior e superior do coeficiente de variação, gruas de liberdade
sagui_mont_tumuc_hn$`Sem termo`$dht$individuals$D
NA
# área de estudo, tamanho da área de estudo, area coberta pelo esforço amostral, esforço amostral em metros, número de detecções, número de transectos (ea), taxa de encontro, coeficiente de variação da taxa de encontro
sagui_mont_tumuc_hr$`Sem termo`$dht$individuals$summary[1:9]
# área de estudo, tamanho da área de estudo, trilhas ou estações amostrais, esforço total em cada trilha, abundância estimada em cada estação amostral, número de detecções em cada estação amostral, área total amostrada
sagui_mont_tumuc_hr$`Sem termo`$dht$individuals$Nhat.by.sample[1:8]
# total, densidade estimada, erro padrão da densidade destimada, coeficiente de variação da densidade destimada, intervalo de confiança inferior e superior do coeficiente de variação, gruas de liberdade
sagui_mont_tumuc_hr$`Sem termo`$dht$individuals$D
Fonte: zoochat.com
cutia_alto_tarau <- transformar_para_distanceR_covariaveis() |>
filter(
Region.Label == "Resex Alto Tarauacá",
sp_name == "Myoprocta pratti"
) |>
drop_na(distance)
cutia_alto_tarau |>
plotar_distribuicao_distancia_interativo()
Warning: Continuous y aesthetic
ℹ did you forget `aes(group = ...)`?
cutia_alto_tarau_hn <- cutia_alto_tarau |>
ajuste_modelos_distance_hn(truncamento = 10)
Fitting half-normal key function
AIC= 549.734
Starting AIC adjustment term selection.
Fitting half-normal key function
AIC= 549.734
Fitting half-normal key function with cosine(2) adjustments
AIC= 551.371
Half-normal key function selected.
Starting AIC adjustment term selection.
Fitting half-normal key function
AIC= 549.734
Fitting half-normal key function with Hermite(4) adjustments
AIC= 551.211
Half-normal key function selected.
cutia_alto_tarau_hr <- cutia_alto_tarau |>
ajuste_modelos_distance_hr(truncamento = 10)
Fitting hazard-rate key function
AIC= 550.863
Starting AIC adjustment term selection.
Fitting hazard-rate key function
AIC= 550.863
Fitting hazard-rate key function with cosine(2) adjustments
AIC= 552.868
Hazard-rate key function selected.
Starting AIC adjustment term selection.
Fitting hazard-rate key function
AIC= 550.863
Fitting hazard-rate key function with simple polynomial(4) adjustments
AIC= 552.868
Hazard-rate key function selected.
cutia_alto_tarau_hn |>
purrr::map(\(.x) plot(.x))
$`Sem termo`
NULL
$Cosseno
NULL
$`Hermite polinomial`
NULL
cutia_alto_tarau_hr |>
purrr::map(\(.x) plot(.x))
$`Sem termo`
NULL
$Cosseno
NULL
$`Polinomial simples`
NULL
summarize_ds_models(
cutia_alto_tarau_hn$`Sem termo`,
cutia_alto_tarau_hn$Cosseno,
cutia_alto_tarau_hn$`Hermite polinomial`,
cutia_alto_tarau_hr$`Sem termo`,
cutia_alto_tarau_hr$Cosseno,
cutia_alto_tarau_hr$`Polinomial simples`
)
cutia_alto_tarau_hn |>
purrr::map(\(.x) gof_ds(model = .x))
$`Sem termo`
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 0.346228 p-value = 0.100677
$Cosseno
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 0.346228 p-value = 0.100677
$`Hermite polinomial`
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 0.346228 p-value = 0.100677
cutia_alto_tarau_hr |>
purrr::map(\(.x) gof_ds(model = .x))
$`Sem termo`
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 0.36764 p-value = 0.0881081
$Cosseno
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 0.36764 p-value = 0.0881081
$`Polinomial simples`
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 0.36764 p-value = 0.0881081
# área de estudo, tamanho da área de estudo, area coberta pelo esforço amostral, esforço amostral em metros, número de detecções, número de transectos (ea), taxa de encontro, coeficiente de variação da taxa de encontro
cutia_alto_tarau_hn$`Sem termo`$dht$individuals$summary[1:9]
# área de estudo, tamanho da área de estudo, trilhas ou estações amostrais, esforço total em cada trilha, abundância estimada em cada estação amostral, número de detecções em cada estação amostral, área total amostrada
cutia_alto_tarau_hn$`Sem termo`$dht$individuals$Nhat.by.sample[1:8]
# total, densidade estimada, erro padrão da densidade destimada, coeficiente de variação da densidade destimada, intervalo de confiança inferior e superior do coeficiente de variação, gruas de liberdade
cutia_alto_tarau_hn$`Sem termo`$dht$individuals$D
NA
# área de estudo, tamanho da área de estudo, area coberta pelo esforço amostral, esforço amostral em metros, número de detecções, número de transectos (ea), taxa de encontro, coeficiente de variação da taxa de encontro
cutia_alto_tarau_hr$`Sem termo`$dht$individuals$summary[1:9]
# área de estudo, tamanho da área de estudo, trilhas ou estações amostrais, esforço total em cada trilha, abundância estimada em cada estação amostral, número de detecções em cada estação amostral, área total amostrada
cutia_alto_tarau_hr$`Sem termo`$dht$individuals$Nhat.by.sample[1:8]
# total, densidade estimada, erro padrão da densidade destimada, coeficiente de variação da densidade destimada, intervalo de confiança inferior e superior do coeficiente de variação, gruas de liberdade
cutia_alto_tarau_hr$`Sem termo`$dht$individuals$D
NA
Fonte: zoochat.com
macaco_jaru <- transformar_para_distanceR_covariaveis() |>
filter(
Region.Label == "Rebio do Jaru",
sp_name == "Lagothrix cana"
) |>
drop_na(distance)
macaco_jaru |>
plotar_distribuicao_distancia_interativo()
Warning: Continuous y aesthetic
ℹ did you forget `aes(group = ...)`?
macaco_jaru_hn <- macaco_jaru |>
ajuste_modelos_distance_hn(truncamento = 20)
Fitting half-normal key function
AIC= 699.333
Starting AIC adjustment term selection.
Fitting half-normal key function
AIC= 699.333
Fitting half-normal key function with cosine(2) adjustments
AIC= 688.909
Fitting half-normal key function with cosine(2,3) adjustments
Warning: Detection function is not strictly monotonic!AIC= 689.303
Half-normal key function with cosine(2) adjustments selected.
Starting AIC adjustment term selection.
Fitting half-normal key function
AIC= 699.333
Fitting half-normal key function with Hermite(4) adjustments
AIC= 701.224
Half-normal key function selected.
macaco_jaru_hr <- macaco_jaru |>
ajuste_modelos_distance_hr(truncamento = 20)
Fitting hazard-rate key function
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).AIC= 629.836
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).Starting AIC adjustment term selection.
Fitting hazard-rate key function
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).AIC= 629.836
Fitting hazard-rate key function with cosine(2) adjustments
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).Warning: Model fitting did not converge. Try different initial values or different model
Model failed to converge.
Hazard-rate key function selected.
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).Starting AIC adjustment term selection.
Fitting hazard-rate key function
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).AIC= 629.836
Fitting hazard-rate key function with simple polynomial(4) adjustments
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).Warning: Model fitting did not converge. Try different initial values or different model
Model failed to converge.
Hazard-rate key function selected.
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).
macaco_jaru_hn |>
purrr::map(\(.x) plot(.x))
$`Sem termo`
NULL
$Cosseno
NULL
$`Hermite polinomial`
NULL
macaco_jaru_hr |>
purrr::map(\(.x) plot(.x))
$`Sem termo`
NULL
$Cosseno
NULL
$`Polinomial simples`
NULL
summarize_ds_models(
macaco_jaru_hn$`Sem termo`,
macaco_jaru_hn$Cosseno,
macaco_jaru_hn$`Hermite polinomial`,
macaco_jaru_hr$`Sem termo`,
macaco_jaru_hr$Cosseno,
macaco_jaru_hr$`Polinomial simples`
)
macaco_jaru_hn |>
purrr::map(\(.x) gof_ds(model = .x))
$`Sem termo`
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 1.56017 p-value = 0.000125973
$Cosseno
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 1.26802 p-value = 0.000587037
$`Hermite polinomial`
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 1.56017 p-value = 0.000125973
macaco_jaru_hr |>
purrr::map(\(.x) gof_ds(model = .x))
$`Sem termo`
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 5.3454 p-value = 0.000122761
$Cosseno
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 5.3454 p-value = 0.000122761
$`Polinomial simples`
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 5.3454 p-value = 0.000122761
# área de estudo, tamanho da área de estudo, area coberta pelo esforço amostral, esforço amostral em metros, número de detecções, número de transectos (ea), taxa de encontro, coeficiente de variação da taxa de encontro
macaco_jaru_hn$`Sem termo`$dht$individuals$summary[1:9]
# área de estudo, tamanho da área de estudo, trilhas ou estações amostrais, esforço total em cada trilha, abundância estimada em cada estação amostral, número de detecções em cada estação amostral, área total amostrada
macaco_jaru_hn$`Sem termo`$dht$individuals$Nhat.by.sample[1:8]
# total, densidade estimada, erro padrão da densidade destimada, coeficiente de variação da densidade destimada, intervalo de confiança inferior e superior do coeficiente de variação, gruas de liberdade
macaco_jaru_hn$`Sem termo`$dht$individuals$D
NA
# área de estudo, tamanho da área de estudo, area coberta pelo esforço amostral, esforço amostral em metros, número de detecções, número de transectos (ea), taxa de encontro, coeficiente de variação da taxa de encontro
macaco_jaru_hr$`Sem termo`$dht$individuals$summary[1:9]
# área de estudo, tamanho da área de estudo, trilhas ou estações amostrais, esforço total em cada trilha, abundância estimada em cada estação amostral, número de detecções em cada estação amostral, área total amostrada
macaco_jaru_hr$`Sem termo`$dht$individuals$Nhat.by.sample[1:8]
# total, densidade estimada, erro padrão da densidade destimada, coeficiente de variação da densidade destimada, intervalo de confiança inferior e superior do coeficiente de variação, gruas de liberdade
macaco_jaru_hr$`Sem termo`$dht$individuals$D
NA
Fonte: zoochat.com
cervo_tap_arap <- transformar_para_distanceR_covariaveis() |>
filter(
Region.Label == "Resex Tapajos-Arapiuns",
sp_name == "Mazama americana"
) |>
drop_na(distance)
cervo_tap_arap |>
plotar_distribuicao_distancia_interativo()
Warning: Continuous y aesthetic
ℹ did you forget `aes(group = ...)`?
cervo_tap_arap_hn <- cervo_tap_arap |>
ajuste_modelos_distance_hn(truncamento = 13)
Fitting half-normal key function
AIC= 532.34
Starting AIC adjustment term selection.
Fitting half-normal key function
AIC= 532.34
Fitting half-normal key function with cosine(2) adjustments
Warning: Detection function is not strictly monotonic!AIC= 520.359
Fitting half-normal key function with cosine(2,3) adjustments
Warning: Detection function is not strictly monotonic!AIC= 700.954
Half-normal key function with cosine(2) adjustments selected.
Warning: Detection function is not strictly monotonic!Starting AIC adjustment term selection.
Fitting half-normal key function
AIC= 532.34
Fitting half-normal key function with Hermite(4) adjustments
AIC= 534.295
Half-normal key function selected.
cervo_tap_arap_hr <- cervo_tap_arap |>
ajuste_modelos_distance_hr(truncamento = 13)
Fitting hazard-rate key function
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).AIC= 451.562
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).Starting AIC adjustment term selection.
Fitting hazard-rate key function
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).AIC= 451.562
Fitting hazard-rate key function with cosine(2) adjustments
Warning: Hessian not computable after method nlminbWarning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).Warning: Model fitting did not converge. Try different initial values or different model
Model failed to converge.
Hazard-rate key function selected.
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).Starting AIC adjustment term selection.
Fitting hazard-rate key function
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).AIC= 451.562
Fitting hazard-rate key function with simple polynomial(4) adjustments
Warning: Hessian not computable after method nlminbWarning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).Warning: Model fitting did not converge. Try different initial values or different model
Model failed to converge.
Hazard-rate key function selected.
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).
cervo_tap_arap_hn |>
purrr::map(\(.x) plot(.x))
$`Sem termo`
NULL
$Cosseno
NULL
$`Hermite polinomial`
NULL
cervo_tap_arap_hr |>
purrr::map(\(.x) plot(.x))
$`Sem termo`
NULL
$Cosseno
NULL
$`Polinomial simples`
NULL
summarize_ds_models(
cervo_tap_arap_hn$`Sem termo`,
cervo_tap_arap_hn$Cosseno,
cervo_tap_arap_hn$`Hermite polinomial`,
cervo_tap_arap_hr$`Sem termo`,
cervo_tap_arap_hr$Cosseno,
cervo_tap_arap_hr$`Polinomial simples`
)
cervo_tap_arap_hn |>
purrr::map(\(.x) gof_ds(model = .x))
$`Sem termo`
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 1.85018 p-value = 2.778e-05
$Cosseno
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 1.45409 p-value = 0.000219796
$`Hermite polinomial`
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 1.85018 p-value = 2.778e-05
cervo_tap_arap_hr |>
purrr::map(\(.x) gof_ds(model = .x))
$`Sem termo`
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 2.99603 p-value = 6.83111e-07
$Cosseno
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 2.99603 p-value = 6.83111e-07
$`Polinomial simples`
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 2.99603 p-value = 6.83111e-07
# área de estudo, tamanho da área de estudo, area coberta pelo esforço amostral, esforço amostral em metros, número de detecções, número de transectos (ea), taxa de encontro, coeficiente de variação da taxa de encontro
cervo_tap_arap_hn$`Sem termo`$dht$individuals$summary[1:9]
# área de estudo, tamanho da área de estudo, trilhas ou estações amostrais, esforço total em cada trilha, abundância estimada em cada estação amostral, número de detecções em cada estação amostral, área total amostrada
cervo_tap_arap_hn$`Sem termo`$dht$individuals$Nhat.by.sample[1:8]
# total, densidade estimada, erro padrão da densidade destimada, coeficiente de variação da densidade destimada, intervalo de confiança inferior e superior do coeficiente de variação, gruas de liberdade
cervo_tap_arap_hn$`Sem termo`$dht$individuals$D
NA
# área de estudo, tamanho da área de estudo, area coberta pelo esforço amostral, esforço amostral em metros, número de detecções, número de transectos (ea), taxa de encontro, coeficiente de variação da taxa de encontro
cervo_tap_arap_hr$`Sem termo`$dht$individuals$summary[1:9]
# área de estudo, tamanho da área de estudo, trilhas ou estações amostrais, esforço total em cada trilha, abundância estimada em cada estação amostral, número de detecções em cada estação amostral, área total amostrada
cervo_tap_arap_hr$`Sem termo`$dht$individuals$Nhat.by.sample[1:8]
# total, densidade estimada, erro padrão da densidade destimada, coeficiente de variação da densidade destimada, intervalo de confiança inferior e superior do coeficiente de variação, gruas de liberdade
cervo_tap_arap_hr$`Sem termo`$dht$individuals$D
NA
Fonte: pinterest.com
inambu_mont_tumuc <- transformar_para_distanceR_covariaveis() |>
filter(
Region.Label == "Parna Montanhas do Tumucumaque",
sp_name == "Tinamus major"
) |>
drop_na(distance)
inambu_mont_tumuc |>
plotar_distribuicao_distancia_interativo()
Warning: Continuous y aesthetic
ℹ did you forget `aes(group = ...)`?
inambu_mont_tumuc_hn <- inambu_mont_tumuc |>
ajuste_modelos_distance_hn(truncamento = 15)
Fitting half-normal key function
AIC= 629.656
Starting AIC adjustment term selection.
Fitting half-normal key function
AIC= 629.656
Fitting half-normal key function with cosine(2) adjustments
AIC= 626.077
Fitting half-normal key function with cosine(2,3) adjustments
AIC= 627.756
Half-normal key function with cosine(2) adjustments selected.
Starting AIC adjustment term selection.
Fitting half-normal key function
AIC= 629.656
Fitting half-normal key function with Hermite(4) adjustments
AIC= 631.591
Half-normal key function selected.
inambu_mont_tumuc_hr <- inambu_mont_tumuc |>
ajuste_modelos_distance_hr(truncamento = 15)
Fitting hazard-rate key function
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).AIC= 578.02
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).Starting AIC adjustment term selection.
Fitting hazard-rate key function
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).AIC= 578.02
Fitting hazard-rate key function with cosine(2) adjustments
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).Warning: Model fitting did not converge. Try different initial values or different model
Model failed to converge.
Hazard-rate key function selected.
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).Starting AIC adjustment term selection.
Fitting hazard-rate key function
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).AIC= 578.02
Fitting hazard-rate key function with simple polynomial(4) adjustments
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).Warning: Model fitting did not converge. Try different initial values or different model
Model failed to converge.
Hazard-rate key function selected.
Warning: Estimated hazard-rate scale parameter close to 0 (on log scale). Possible problem in data (e.g., spike near zero distance).
inambu_mont_tumuc_hn |>
purrr::map(\(.x) plot(.x))
$`Sem termo`
NULL
$Cosseno
NULL
$`Hermite polinomial`
NULL
inambu_mont_tumuc_hr |>
purrr::map(\(.x) plot(.x))
$`Sem termo`
NULL
$Cosseno
NULL
$`Polinomial simples`
NULL
summarize_ds_models(
inambu_mont_tumuc_hn$`Sem termo`,
inambu_mont_tumuc_hn$Cosseno,
inambu_mont_tumuc_hn$`Hermite polinomial`,
inambu_mont_tumuc_hr$`Sem termo`,
inambu_mont_tumuc_hr$Cosseno,
inambu_mont_tumuc_hr$`Polinomial simples`
)
inambu_mont_tumuc_hn |>
purrr::map(\(.x) gof_ds(model = .x))
$`Sem termo`
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 0.955575 p-value = 0.00312854
$Cosseno
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 0.820002 p-value = 0.00655342
$`Hermite polinomial`
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 0.955575 p-value = 0.00312854
inambu_mont_tumuc_hr |>
purrr::map(\(.x) gof_ds(model = .x))
$`Sem termo`
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 4.92203 p-value = 6.83418e-05
$Cosseno
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 4.92203 p-value = 6.83418e-05
$`Polinomial simples`
Goodness of fit results for ddf object
Distance sampling Cramer-von Mises test (unweighted)
Test statistic = 4.92203 p-value = 6.83418e-05
# área de estudo, tamanho da área de estudo, area coberta pelo esforço amostral, esforço amostral em metros, número de detecções, número de transectos (ea), taxa de encontro, coeficiente de variação da taxa de encontro
inambu_mont_tumuc_hn$`Sem termo`$dht$individuals$summary[1:9]
# área de estudo, tamanho da área de estudo, trilhas ou estações amostrais, esforço total em cada trilha, abundância estimada em cada estação amostral, número de detecções em cada estação amostral, área total amostrada
inambu_mont_tumuc_hn$`Sem termo`$dht$individuals$Nhat.by.sample[1:8]
# total, densidade estimada, erro padrão da densidade destimada, coeficiente de variação da densidade destimada, intervalo de confiança inferior e superior do coeficiente de variação, gruas de liberdade
inambu_mont_tumuc_hn$`Sem termo`$dht$individuals$D
NA
# área de estudo, tamanho da área de estudo, area coberta pelo esforço amostral, esforço amostral em metros, número de detecções, número de transectos (ea), taxa de encontro, coeficiente de variação da taxa de encontro
inambu_mont_tumuc_hr$`Sem termo`$dht$individuals$summary[1:9]
# área de estudo, tamanho da área de estudo, trilhas ou estações amostrais, esforço total em cada trilha, abundância estimada em cada estação amostral, número de detecções em cada estação amostral, área total amostrada
inambu_mont_tumuc_hr$`Sem termo`$dht$individuals$Nhat.by.sample[1:8]
# total, densidade estimada, erro padrão da densidade destimada, coeficiente de variação da densidade destimada, intervalo de confiança inferior e superior do coeficiente de variação, gruas de liberdade
inambu_mont_tumuc_hr$`Sem termo`$dht$individuals$D
NA